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ABSTRACT 


A three-level quasi-geostrophic atmospheric model on a 
mid-latitude 8-plane with a single wave in the x-direction 
and a thick lower layer was developed. This model was used 
Lo determine tbe disturbance structure in y for linear inter- 
actions with a constant mean wind and nonlinear interactions 
with a varying mean wind. Simple heating to maintain the 
I-neundgcanóosbpeseson were included in the nonlinear case. 

The linear experiment produced exponential growth rates 
in the disturbance field that were independent of the lower 
layer thickness. These waves were baroclinically unstable. 
ic nean wid tela showed large modifications for a Rossby 


deformation length of 75 KM at 104 Jupiter days. 
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IESZEENTRODDETTON 


The meridional belts (dark bands) and zones (light bands) 
of Jupiter have been recorded since this planet was first 
EB ervegd-chrouehza telescope. The features found in these 
regions are semi-permanent and have a variety of rotation 
pcs | Peek (10956) Assembled this information to show the 
variation of zonal velocity with Jovian latitude. Evidence 
was found by Chapman (1969) for anticyclonic shear to be 
associated with the zones and for cyclonic shear to be 
Eco a Sa but the correlation was imperfect. 

Ingersoll and Cuzzi (1969) have investigated the baro- 
tropic instability of che observed motions at mid-latitudes 
using data from Peek (1958). Geostrophic balance and 
symmetric temperature difference between the belts and zones 
neveras media tne Validity of the assumptions was verified 
by the excellent agreement between predicted and observed 
A T it was also shown that the relative 
uem d gradient U" approaches, but does not exceed the 
planetary vorticity gradient 8. This is the necessary 
condition for Bamerrophie instability. Thus barotrophic 
instability probably limits either the belt wave number 
Caem Compile of the Zonal flow. It was further shown 
that the zones are warmer cnans the Delts. 

Stone (1967) applied the baroclinic stability theory 


Redy namics of the Jovian atmosphere. He tested the 





hypothesis that the zonal motions in Jupiter's atmosphere 
are thermal winds and that the latitudinal cloud bands are 
caused by baroclinic instabilities under non-geostrophic 
conditions. Stone (1966a) showed that at least three 
Ate rent Kinds Of baroclinic instability can occur when 
de/dZ » 0. The kind of instability that dominates depends 


upon the value of the non-dimensional Richardson number 


_ & 39 ,,8U.* 
ri = & 28/2)", (1.1) 


Da ccesrussbheNaceeleration of gravity, 6 às the potential 


temperature, = -— is a representative vertical stability, 
MR tic vertical shear of the wind. The 


dZ 


aras operas a characteristic ratio of work done 
against gravitational stability to energy transferred from 
Mel cOncuUropuUrene MOtion. He assumed that the observed 
(Comoe we meer are due to Symmetrical instabilities with 

Io cian variations. This type of instability domin- 
ates in the range 0.25 < Ri < 0.95. Allowance for a deep 
atmosphere was made by not using the Boussinesq approximation 
for compressibility. It was shown that in a deep atmosphere 
symmetric instabilities form in the specified range of Ri 

and that these instabilities may extend over many scale 
heights. Thus the zonal winds and the latitudinal cloud 
bands may be explained qualitatively by assuming the presence 
of north-sourth temperature gradients caused by differential 


solar heating. These gradients will drive thermal winds in 





DM rection, and under appropriate conditions, 
Eu amstabälities can grow, generating latitudinal 
etoudy bands. 

Gierasch, Ingersoll, and Williams (1973) demonstrated 
Bias a cloudy planetary atmosphere at rest is unstable to 
disturbances of large horizontal scale. The energy source 
Bscgsemsrability is the change in radiative heat flux 
Hechos svercical displacement near the emitting 
level. The important feature of the model was that the net 
hcc ave Qs iproportional to the vertical displacement 
DN neselosudeEnPom rts equilibrium position.  Quasi-steady 
geostrophic MAR So short period oscillations, representing 
inertia-gravity waves may become unstable in the presence 
of radiative heating. Gierasch and Goady (1969) found that 
Pc a mo conScanits of Jupiter are longer than the 
dynamical time constants in the earth's atmosphere; there- 
uo 1. Oph e Mode was applied, since the instability 
can be quasi-geostrophic. Gierasch, Ingersoll, and Williams 
(1973) derived a value of 19,000 KM for the y-scale wavelength 
QUNM latitude Lus corresponds to the limiting amplitude 
for an axisymmetric geostrophic instability where the wave- 
ene x-direction Ky cu UIS y- length 15S essen- 
tially equal to the observed belt-zone pair at this latitude. 
They also computed the value of the radiative time constant 


X -1 


, For an isothermal atmosphere E was approximately 


equal Lc 160 years. It was also observed that even if the 


lapse rate differed from the adiabatic lapse rate by only 





12, the radiative time constant would be about 1.6 years. 
ME cceucoculveSSarceconsrstenb with the observed lifetimes 
RaEctHospheryo features on Jupiter. 

Gierasch (1973) studied the radiative instability of 
a cloud deck on an equatorial B-plane as discussed by 
Gierasch, Ingersoll, and Williams (1973). The model was 
based on the assumption that a precipitating condensable 
atmospheric constituent can be concentrated in the atmos- 
RR re by vertical motion. If the condensing constituent is 
radiatively important, it can cause radiative warming at 
Migem ocartions where dt iS concentrated. As a result, the 
heavtine is related Co vertical motion and reinforcement is 
possible. They considered two values of the Rossby 


deformation length 


Beyer E (1.2) 


where g' is reduced gravity, f is the Coriolis parameter, 
EMEN sa characteristic height. Static stability esti- 
mates by Lewis (1969) and Stone (1972) gave characteristic 
lengths of 785 KM and 78.5 KM, respectively. It was con- 
cluded that the radiative instability mechanism driven by 

the dynamical concentration of a condensable constituent, 
produces instabilities with the wavelength and zonal symmetry 
of Jupiter's bands. Nonsymmetric instabilities with a 


shorter wavelength also exist. Finally, the structure of 
\ 





the modes is in agreement with Ingersoll and Cuzzi's (1969) 
conclusion that cloudy zones are warmer than belts. 

ce = level quasi-seostrophie model incorporsting a 
deep lower layer, as suggested by Gierasch (1973), was 
developed based on the Charney and Phillips (1953) multi- 
level model. Tne model was modified by implementing Lorenz's 
(1960) technique of replacing the dependent variable with 
Fourier Ao cents. 

The objective of this study was to define the linear 
disturbance structure in y and to see if the belt-zone 
wind growth rate due to quasi-geostrophic baroclinic 
instability decreases with Re lower layer depth. Non- 
linear interactions with the belt-zone flow pattern were 
eRenzstndted aáncorporatine simple heating, to maintain 
IT TG and friction. These objectives were 
evaluated using two values for the Rossby deformation 


length. 





II. THE MODEL AND FORECAST EQUATIONS 


A. DESCRIPTION 

The atmospheric model for these experiments was based 
on the quasi-geostrophic equations. The model was constructed 
on a B-plane with a central latitude of 45 degrees. Since 
the large-scale motions in the extra-tropical latitudes of 
Jupiter are quasi-geostrophic, the domain of the model was 
restricted to the middle latitudes. Phillips (1963) dis- 
cussed geostrophic motions with approximately uniform 
Perential vorticity. These motions have characteristic 
wavelengths of about 4,000 KM on earth, or the order of 
the Rossby deformation length. He used an expansion in 
the Rossby number, Rz U7 tie waeresU is A characteristic 
Be locıty, f is the Coriolis parameter, and L a character- 
istic length. This expansion showed that the f-plane 
approximation is justified for geostrophic motions of this 
type. Gierasch, Ingersoll, and Williams (1973) have shown 
that large-scale bands have a characteristic scale of 15,000 
KM in the atmosphere of Jupiter. It was further shown that 
the Rossby number Ro is much less than one. If the rota- 
Coe e ciently fast such that p 1S Ol Order Oneto 
less, Coriolis forces become important and must be considered 
im che equations of motion. The additional conditions of 
R,/a, where a is the planetary radius, being small, and the 


large-scale waves being greater than or equal to Rr, are 


10 





E rTSed. vhepbrefore, the geostrophic assumption uc. 
NHucEScophystreat constants for earth and Jupiter are summa- 
red in the appendix. 

eree- level model was developed so that heating and 
ERON could be included and still satisfy the lower 
Binary condition. Since the interior of the planet is 
Mco, a deep lower layer is used to remove the two 
Peers ayers from the lower boundary. The solution function 
ENDE US derivatives were forced to vanish at the boundaries. 
The Jovian atmosphere was divided into six layers numbered 
1 through 6 as shown in Figure 1. The upper layers, 0-4, 
have thickness AP/2; and the lower layers, 5 and 6, have 
erreiness P/2 here P>>AP in this application. The thick- 
Mess between levels 3 and 5 is defined to be mn = (AP+P)/2. 


in ceneral these thickness relationships are arbitrary. 


Pe METROD OR SOLUTION 


The quasi-geostrophic vorticity equation on a -plane is 


Em. du. po? 
TE tV V(EtB.y) - T. 5 pU (2 1) 
Waere g = = V6 , V= = kxVo , (22) 
O O 


and the thermodynamic equation is 


d d 


+ t VV + Wo = E (2.3) 
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(397 90 
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FIGURE 1 - The Three-Level Model 
Variablesare evaluated at the layers shown. 
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where u = -- and o is defined as the static stability. The 
. A A us : o 

weyn 1S Che radiative time constant and SP is an 

Dan temperature which is dependent upon radiation, 


convection, and other small scale parameters. Equation (2.1) 


Ras applied at leveis 1, 3, and 5 which yielded 








304 “ao 2 

RE + NA E L. E = BV Ch > (2.4) 
dL 0, —U) 

955 | e 2 6 
or ICM uev E ls Pe ze) LE (2.6) 


Equation (2.3) applied to levels 2 and 4 yielded 


X 
ob d 
TONER 


93791 
Poss 


Ta 
Pa-Py 








d 
em 


Ale 


) + V.-Vf O 
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CTO 


2 ÓN 


Q2 


d C ^ 2:019 
— ( ) + Vy V É s 940] 


cu E 





96" 96 
S - GP 
(2.8) 


Q2 
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Equations (2.7) and (2.8) were solved for o, and wys 
respectively, and substituted with Equation (2.2) into 
Equations (2.4), (2.5), and (2.6). These operations yielded 


the potential vorticity equations 


15 





2 2 z .y(y? 
ES LV $4 *k5(65704 2] B “Vy V(V $, *f B y) = K5V5* V(65764) 


Q2 


4 


K x 
+ BV 6) + —365-(05-6,2] - €, , (2.9) 


[V^ 6, *k4(65-64)-k5(0570,)1 9 -V,- V(V^ o, f B y) 


a] a 
er 


= KT: VL dd) + K5V5* VC 62764) 


K 
* BV 6. * 2201-65-04) 


És x u 
Ex dos = (dad) = Co , (2. 


2 - d 
= [y $5-Kg ($57 64)] = zu UY $e tf Boy) d KgVy* V(05765) 


y E u 
* BV $5 - —364-(05-9,)1 7» €, , (2. 











qe | 
where Ka = : qo 
2 2 
3 AP 05 
TAS 
Ky = = 6K. 3 G2. 
APTO y, 
e 
and Epis = Ek) . (2. 
ROST 


14 


10) 


IL) 


12) 


13) 


14) 





A set of three linearly independent eigenvector equa- 
tions were formed to solve the equation set (2.9), (2.10), 
ES uUsPter "a method by Charney and Phillips (1953). 


Er 


2 3 (2515) 


then the general form of the potential vorticity equation 


ES 
2 90 o0 — 
V 3t T A “a JE x C 3 (2.06) 
Cy 
where o> Co , 
3 
cs che appropriate matrix. If 
iL 
u =| > 2) 
Y 


X 


is an eigenvector of A and A is its corresponding eigenvalue, 


then 
Au = Au : 21.5) 
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and 


-] JL U 
A = Į -(1+8) E (2. 19) 
U 6 -€ 
Now if 
o = H y + Hopo + Ham : (2220) 
and C = du + dou, + dom, ; (2.210) 
Re Substituted into Equation (2.16) and Equation (2.18) is 


PE oe e iimne che eigenvalues, then the coefficients 


each eigenvector can be equated to give 


ie 
at 


Mere i = 1,2,3. 


The following eigenvalues as specified by Equation (2.18) 


were obtained, 


[(v 43, k4)H, ] - d, , (2. 


of 


22) 


p «8 
ELE > CATA > [(e45)545 (1-6) ]* (2.23) 
Ay mre l (2464€) E 5 [Cers)2+1(1-€)]% 
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Ine eigenvector equation set was formed from the poten- 
ELM eey equations (2.9), (2.10), and (2.11) and 


Equation (2.22). These equations are 


3 2 me. 2 2 26 
mU [V H, ecu [8 YKxVH, +V, V $1 *V4V DE a5 J 


at 
4 By* 
VH, + Q : (o TNI 
9 2 E | 2 2 2 
zg LV *14k5)E,] - -V (BoYkxVH5*V, V $4 *85V.V datyoVEV és 
= Às K o LV, (937 $1) "YaYy(d5- $3421) iE By Hs T Q 3 
(orais) 
9 E 2 
SE ES 32 Ho ] VB ¿YKXVH +7 vê $1 *84V.V "paT a VEN óc 
lj 
- Agk«5lV5 (ó.-6,)-Y5V4(0570.)]) + BV H. +Q , 
(2.26) 
E a x x 
where Q = - LüG-8)(805-8305) + (38-Y€) (39-96 1)] ; (2:205 
An assumed solution of the form 
H = H(y,t) + H¿cos kx + H, sin kx , (2525 


where H is the zonal mean solution amplitude, Ha Mc iviee 
Fourier cosine eddy amplitude, and Ho Poem Outer sine 
eddy amplitude, was introduced into Equations (2.24), (2.25), 
and (2.26) as suggested by Lorenz (1960). The equations 


were quasi-linearized by eliminating all terms of higher 
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wave number than k. This technique permits full interaction 
oí the disturbance in x with the mean flow, but neglects 
interaction with the other waves in x. 

BOuacions (2.24), (2.25), and (2.26) were separated 
No components independent of x, coefficients of cos kx, 
EmenNcociincients oi sin kx. The resulting equation set 
(2.29) through (2.43) constitute the prediction equations 
for this model. The primed quantities represent the partial 
ucpo)vative with respect to y. 


The terms independent of x in (2.24) give: 


3 tt e k 11! ttt 11 ] "n t 
n° BF, | rcbis 91591577 (01 5915703 093,) 

TL Haa "has Daa lL bas Taa "has has 

n 6 111 tt! 6 " t 11 11 

€ Os ex E ee E 

+L N L Q ; (2529) 

where S ON , ELO 
"" E * — * 

and SCC [(1-8)(865-3045)*(88-ve)(30,-30,)1 . (2238) 
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The cos kx terms in (2.24) give: 


ERU 5. -. kn eos GI L k^$, 6, .-9, 6.) 
ot 3y* c Ji ls E 
E (ONE SNC EN 
3 SO as 
> EET: óc DIt +g 
2 5 ROS 9 QUE 
(232) 
where ds B(H, -2K2H,+k E) . (2033) 
k5 x Xx 
and Da « L(-8)(265,-3645,)*(68-y£) (30-20, ,)] 230 


thesin kx terms in (2.23) give: 


2 T 


d d p Z k _ım 2! E 
ot E = - FB of ol B Coy ES E 
E, 9393030) 
ô 111 F. 4 
es (65 ur “E D Ro $. 0.) T Q o 
(2.35) 
where D B(H. -2k^H 4k H ) (2330) 
S S S S > 
> x x 
and Q, = — L(1-8)(365-36,.)*(68-ye) (30,-30,,)] .(2.37 
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The terms independent of x in (2.25) give: 


2 
3 3 > E EN m o RE E " t E " t 
dt ET = AK, JH m ¿fo e ee tee DE o) 
+B2L (92.337 3593 027 (0329307030935) 
A. ese! te tete) 
t t ( t 
apa, Laha" has Haa TL Das h Tac 
AT UA AE 
L N L Q l) (223.0) 
The cos kx terms in (2.25) give: 
2 
Q Q 2 n k Lg po a " 
EV EH c a RR E Dr) 


(BALE, ba tk Faba Tala 


tt Da 11 
TER Pr PK 95 RUD 


* Agk5[(9165,-050, .29Y5(0670,.-0.07.0]) + Fo + Go 
| (539) 
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The sin kx terms in (2.25) give: 


Q»| e» 
ct 


2 


dy 


d 
NE -k 


PAARE, = GE TELHA HT Td) 
E 


+85 (04 Oct Td ao Tad ac 


Bast ETE 


raU, dt Fett) 


SR $564) t Y4 ($564, -9 $; we + Fo + Q. 
(2.40) 


The terms independent of x in (2.26) give: 


2 
(2 + 
ay 


a ie 


3 2) ar, E ts ie) Cis ter fic T 


O “(3 50307F 30036) ] 
HE (5 0%567¢ 56950) 7 (05,65, 765, 65,21 
Eo (939137035010 =(07 09357975030) 
- Y (6 663.765,05, +13 (0300557035050) 1) 


+ F+Q , (2.41) 
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cos kx terms in (2.26) give: 


2 
Q 2 k ett Je d 
E 38518 um (85 f SH. y (dy EU 91915 91915) 


f! 


Io) 
ct 


rn 2t eu 2 pnt om uL. oM 
Ha-pa fk. 0503,70403, *Y3 (65 Oss HE $5 055-95 955) 


+AzkoL (97 63570307 5)+13 (9593570305521) + Fa D da > 


(2.42) 
The sin kx terms in (2.26) give: 
Q (32 T Ts k M m an On 
Re scr Q N E H OT 


— tff 2 t tt 
EB, 002 07.78 9202079302.) 
e C AME 
ER C Bench GC 


W e a UOS 
(2.13) 


A method developed by Richtmeyer (1957) was used to 


o , at each time step. The equa- 


tions were integrated in time by a marching process. All 


solve for the tendency, 


finite differences were centered in time or space. A 


forward time step was used to start the marching process. 


A finite difference scheme developed by Matsuno (1966) was 
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EI Very 50 time steps. This is a two step iteration to 
simulate the backward difference method. This scheme was 
ed an this situation to suppress the computational 
Dede that would arise from using the centered time differ- 


encine exclusively. A time step of 30 minutes was used. 
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SEES UIS 


Two experiments were conducted under varied conditions. 
erste experiment was designed to yield the growth rate 
of the most unstable wave in the disturbance fields for 
Rossby deofrmation lengths of 75 KM and 750 KM. This range 
of values for the Rossby deformation length was suggested 
by Gierasch (1973). The upper two levels each had a pressure 
we erential of 50 CB, while thickness of the lowest level 
= med values of 50, 500, and 5000 CB, respectively. The 
y-wavelength of the mean flow was 15,000 KM for all cases. 
The x-wavelength varied from 1,000 KM to 15,000 KM for 
= SO KM and from 300 KM to 3000 KM for R 


L L 


K, as defined by Equation (2.12) was scaled to be 


= 75 KM. 


K, = —. (315) 
he static stability at level 2, 05, Was ains aaron 
Equation (3.1) aná the static stability at level 4, O y» 
was found by the relation 

E 2 

On = 0, E 3 (ne) 

where P, and Py are the absolute pressures at levels 2 and 


4, respectively. This ratio was a constant equal to 0.25 


for all experiments. The mean wind Us vas initialized as 


24 





Ze museidal function of y with an amplitude of 20.5 
meters per second in level 1. This value was obtained 
from Ur = 2 Nose nietas Value Lor barotrophic 
stability as suggested by Ingersoll and Cuzzi (1969). The 
mean wind was zero in levels 3 and 5. The mean field was 
pouxcounsconecswhylehomakes Lhe equations linear, throughout 
the first experiment. The disturbance field cos kx contained 
an initial sinusoidal variation with the same y-scale as 
the mean flow and an amplitude an order of magnitude smaller 
than the mean field. The initial sin kx field was entirely 
described by the cos kx terms. There were no heating or 
Mate tion terms included. 

The forecasts were carried forward until the growth 
rates became exponential. Growth rates of the disturbance 


were calculated by 


1 939A = ET SA E 
nã NS Mur, ers 


where M is the exponential growth rate, A is the amplitude 

of vis disturbance, and Us and t1 are the times at which 

the disturbance amplitudes were A, and Br > respectively. 

The maximum growth rate for Ri =A OFM occurred at an 
x-wavelength of 6,000 KM for all lower layer thicknesses 
(Figure 2). For lower layer thicknesses of 50 and 500 CB 

and Rr = 75 KM the maximum growth rate was at an x-wavelength 
of 600 KM. The 5000 CB case showed maximum growth rate at 


Fer wexeleneth of 500 KM (Figure 3). 
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GROWTH RATE (105 SEC'!] 
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3 5 7 9 11 13 15 
X-WAVELENGTH (105 METERS] 


FIGURE 2 - Growth Rate of Y-Disturbance Wave for 
Rp=750 KM 
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GROWTH RATE (107^ sEC ^!) 





0 6 12 18 24 30 
X-WAVELENGTH (10? METERS) 


FIGURE 3 - Growth Rate of Y-Disturbance Wave for 


R, = (5 KM 


2 





The phase and amplitude relationships of the disturbance 


field were determined by 


a cos kx + >. Sin kx = © cos(kx — 6), 
sea) 
where o and De are he ecopetentials En the cos kx and 
sin kx fields, respectively. Equation (3.4) was expanded 
to yield the following expressions for the amplitude C 


and the phase angle 9, 


© 
6 = arctan (>) l E 
C 
E Pe s 
C Sos $- ^ sin $0 c i2 


The amplitude and phase relationships for q and 


1? $2» 
$5 were computed at the maximum growth rate for each case. 
The results are shown in Figures ! through 15. A11 ampli- 


tude profiles showed greater relative growth in the northern 


wave; however, the effect was more pronounced at Rr = 750 


KM. This was due to the Peesemee wor sc hezplanebary vorticity 


gradient B. The relative amplitudes of Hu $3» and v 


were the same order of magnitude for the lower level thick- 


ness P = 50 CB (Figures 4 and 7). The amplitude of Po 


was an order of magnitude smaller than $4 or Pa for 
P= 500 CB (Figures 5 and 8). At P = 5000 CB (Figures 6 
and 9) the amplitude of Ne was negligible compared to $4 


or $.. These characteristics were observed at both values 


3 
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FIGURE 4 - Relative Amplitude of Du $5; and $5 
as a Function of Y for 4Ps50 CB, 
R, - 75 KM, and x-wavelength of 600 KM. 
Y-Axis has been shifted to the north. 
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RELATIVE AMPLITUDE 





0 3 6 9 12 15 
-Y (10? KM) 


FIGURE 5 - Relative Amplitude of Py > E and js 
ass ae nun ons AAN SO CE, 
P = 500CB, R = 75 KM, and x-wavelength 
of 600 Km. 
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RELATIVE AMPLITUDE 


Qi Qı 
À Ps A D3 


3 6 9 12 15 
- Y (107 KM) 


FIGURE 6 - Relative Amplitude of $4; $3; and $5 


asar culona 50 CB, 
P= 5000 CB; Rr = 75 KM, anà 


x-wavelength of 600 KM. 
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FIGURE 7 - Relative Amplitude of by > Ha: and $5 
aS com Une eLousOney wien (AP = 50 CB, 
P = 50 CB, Rr = 750 KM, and 
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of the Rossby deformation length. The phase angles for 
each level showed that all of the waves were barotropically 
Stable, since the phase lag in y was in the same sense as 
the horizontal wind shear. Also, the vertical tilt was in 
Ene opposite sense to the vertical shear as is to be ex- 
pred for baroclinic instability. Note that the y-scale 
of the disturbance is much smaller for R, = 75 KM compared 


L 


to Rj = 750 KM. This result was predicted by Stone (1969) 


who showed that in a slowly varying wind field the y-scale 
of baroclinic waves would be of the order of the Rossby 
deformation length. It was observed that the amplitude 
sche disturbance fields for Br = (5 KM was contained in 
two narrow regions centered about 2600 KM and 7500 KM. 


Friction and heating were added to the model in the 


second experiment to study the nonlinear effects. The 


39* 
dD 


the initial value of zs . This has the effect of main- 


equilibrium temperature profile was obtained from 
taining the mean wind, since the heating will tend to 
restore the mean temperature field to its initial conditions 
when the disturbance amplitude is sufficient to cause a 
departure. Heating was present in the mean field only. 

The radiative time constant t was 10 Jupiter days. The 
friction term was present in all fields. The friction 
coefficient B was equal to 1.0 x 10° Snide. 0 x 10° ioe 
Rossby deformation lengths of 750 KM and 75 KM, respectively. 
Since the grid mesh distance was 150 KM and the x-wavelength 


was 600 KM, the grid may have difficulty resolving the wave; 
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phesefore, a larger value of the friction coefficient was 
needed to maintain stability in the model. The mean field 
was allowed full interaction with the disturbance fields. 
ime cases corresponding to the Rossby der rormaction lengths 
Were performed with a lower layer thickness of 500 CB. 

Fill otner initial conditions were as described in the linear 
experiment. 

The model was integrated for 104 Jupiter days with a 
Rossby deformation length of 75 KM. At this time the model 
Became unstable. The corresponding mean zonal wind profile 
and disturbance amplitudes and phase relationships are 
shown in Figures 16 through 18. The maximum wind speeds 
were 25.9 meters per second for the northern easterlies and 
26.3 meters per second for the southern westerlies. Zonal 
velocities also developed in layers 3 and 5. The relative 
pua associated with the maximum wind is the result of 
paecampiitude of the disturbance fields. These amplitudes 
Enc Similar to the linear case except the peaks are now 
Bemodal anc the profile has a broader base, indicating the 
Gisturbance energy has been distributed over a larger 
region. The phase angles showed baroclinic instability 
in the vicinity of the maximum mean zonal wind in the 
Porch ana the south, since the sense of the vertical wind 
shear was opposite the tilt of the phase. 

Because the 750 KM Rossby deformation length case has 
a much slower disturbance growth rate, the experiment was 


run for 536 Jupiter days. The mean zonal wind profile is 


He 
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12 


15 





-20 ZU! 0 10 20 
ZONAL WIND SPEED (METERS PER SECOND} 


FIGURE 16 - Mean Zonal Wind for AP = 50 CB, 


P = 500 CB, R, = 75 KM, and 
x-wavelength of 600 KM for 
Jupiter Day 104. 
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shown in Figure 22. Relatively small changes have occurred 
compared with the other nonlinear case. The mean wind 

speed is 20.7 meters per second in level 1. Small amplitude 
winds have developed in levels 3 and 5; however, no appre- 
ciable wind is observed in level 5 in are south E hlvure 21 
shows the amplitude of the most unstable wave in the distur- 
ce field as a function of time. The growth rate is 


exponential with value 


= SA m 0.21 x 107° sec” 


Q? 


thus the field has not reached quasi-equilibrium. The 
disturbance amplitudes shown in Figure 19 indicate the 
effects of the planetary vorticity gradient B have been 
reduced. The northern disturbance has a relatively greater 
amplitude in all three levels. The relative amplitude of 
us is an order of magnitude smaller as observed in the 
iamear case. 

One inconsistency was present in this experiment which 
mas not been resolved. The growth rates observed in the 
linear mode were an order of magnitude larger than the 
nonlinear mode. This was not a physical result, but an 


error In the program. 
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IV. CONCLUSIONS 


A three-level quasi-geostrophie atmospheric model has 
been examined on a B-plane with a single wave in the x- 
direction and a deep lower layer of varying thickness. The 
disturbance structure in y was investigated to determine 
the linear interactions with a constant mean wind and for 
nonlinear interactions with a varying mean wind. Simple 
heating to maintain the mean wind and friction were included 
in the nonlinear case. Two values of the Rossby deformation 
length were included in the experiments. 

The linear experiment showed that the growth rates in 
the disturbance field are largely independent of the lower 
layer depth (Figure 2). These growth rates are a function 
nu he Rossby deformation length and the x-wavelength. The 
waves were barotropically damped in this case, since the 
Relative vorticity gradient was less than or equal to the 
Miranevary vorticity gradient. This was obvious since the 
phase tilt was in the same sense as the horizontal wind 
shear. 

The heating function was not important in the nonlinear 
experiment, since the model never reached a long term steady 
state. The disturbance fields were barotropically damped. 
In the cases studied the zonal disturbances are baroclinically 
us-e2ble- rf the actual wind structure on Jupiter is baro- 


clinically stable as suggested by Ingersoll and Cuzzi (1969), 


51 





then this could be accomplished by distributing the mean 
emamge through great depth. This would give a very small 
Em ou Whicn could be stabilized by B. The growth rates were 
ERDectea Lo decrease with increasing lower layer depths, 
Bern fact they did not. Thus it may be that the lower 
Maver was not properly formulated by this model. 

cher investigation of the disturbance growth rates 
EN red TO resolve the discrepancy between the linear 
psuenonjineapbr cases. It is also recommended that the grid 
mesh be reduced to more realistically model the 600 KM 
x-wavelength. This should allow use of a smaller friction 


coefficient. 
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APPENDIX 


GEOPHYSICAL AND ASTRONOMICAL CONSTANTS OF JUPITER AND EARTH 


PARAMETER Ju Iren EARTH 
Perihelion 

Distance, A.U. 8 4.951 0.9833 
EN U.-1.496 x 10 KM) 

Amphelion 5.455 OLOT 
Distance, A.U. 

Siderial 10,759.20 365.256 
Period, (day) 

Mielination of 3904! QOO 
Brit tO ecliptic 

Mass, (GM) 1.899 x 10% 5.976 x 10*! 
Density, (GM/CM?) ME DEBIT 
Mean Radius, (KM) 69,828 EE ST 
Equatorial 

Gravitational > 2587 979 
Acceleration, (CM/SEC ) 

Albedo Qr T se 0.36 
Orbital 0.0483 02.017927 
Eccentricity 

Mean Temperature 130° K 330° K 
Solar UD E 
Constant, (Ly,MIN) 

SEO ) 17.6 x 107? 7.3 x 107? 
Observed Mean 50 15 
Wind, (M/SEC) 

Rossby Number, R, 0.0041 0.032 
Rossby = Mb, 19.5 1,000 


deformation Length, (KM) 


D2 


PARAMETER 


Coriolis Parameter, (45°) 


Planetary Vorticity 


Gradient, (45º), (SEC 


mt) 


R, /a 


JURTTER 


2.5 x 10. 


3.692907 


0205 
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H 


12 


SEC 


-1 


EARTH 


OO 


EOS. 


0.63 


SEC 
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